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The Time Invariance Principle, Ecological (Non)Chaos, and 
A Fundamental Pitfall of Discrete Modeling 

Bo Den^ 

Abstract: This paper is to show that most discrete models used for population dynamics in 

ecology are inherently pathological that their predications cannot be independently verified 

' by experiments because they violate a fundamental principle of physics. The result is used 
O ' 

O I to tackle an on-going controversy regarding ecological chaos. Another implication of the 
(N 
H . 

J2 ■ As a result it suggests that researches based on discrete modeling must be closely scrutinized 
and the teaching of calculus and differential equations must be emphasized for students of 

■ biology. 

^ , 1. Introduction. No models in ecology are better known than the Logistic Map, or have played 

Q ■ a greater role in the development of the chaos theory ( [|29l[T9ll20l [3ll27l). Surprisingly, however, 

'F^ ■ there is not a greater controversy than what was generated by the model's prediction that one- 

qh| species populations are inherently chaotic. 

The key prediction of the Logistic Map, = Q{xn, r) := rx„(l — says that increase 

^ ! in the intrinsic reproduction rate r leads to chaotic dynamics for the population x„. Contradicting 

22 . evidence existed before the chaos theory was popularized in ecology. For example, in a 1971 

O ■ study of an aquatic system, McAllister and LeBrasseur ( OlH ) showed that enrichment led to stable 

(N ■ 

(3 ■ equilibrium. Extensive search for field chaos came up equally empty-handed. For example, well- 



established geographic patterns on microtine species ( IfTSl [T6l ) showed that ecological systems 
^ I tend to stabilize down the north-to- south latitude gradient, correlating well with the ultimate energy 

r\ , abundance from the Sun towards the equator. The most comprehensive hunt for ecological chaos 

CT. was down by EUner and Turchin ( ifTSl ). They used 3 different Lyapunove exponent estimators on 

^ ■ a large collection of empirical data. Out of their 21 field data sets, not a single set scored a positive 

^ ' Lyapunove exponent by two of the 3 estimators. Out of their 20 lab data sets, only two scored a 

o3 I positive Lyapunove exponent by two estimators. The inescapable conclusion is overwhelming — 

ecological chaos is not to be expected in the wild. (Although laboratory chaos is possible with 
stringent setups, such systems are never simple. In fact, the dimension required is 3 or higher, c.f. 

Em.) 

The glaring irreconcilability between the theory and reality can only lead to one logical con- 
clusion: the theory is wrong. Otherwise, ecology would be the definitive branch of science that 
logic imperative would have failed. The purpose of this paper is to make a case that the Logistic 
Map and most discrete models used in ecology and life sciences cannot be models for any physical 
process, population dynamics in particular, and their predictions cannot be independently verified 
by experiments. 
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2. Time Invariance Principle. This conclusion rests on a fundamental principle of physics held 
since the time of Copernicus in the 15th century that a physical law should be the same anywhere 
and anytime in the universe. In other words, a law must take the same mathematical form, derivable 
from experiments carried out at independently chosen times and spaces. As a result, the mathe- 
matical formulation of a law must be endowed with such time invariance property. Taken to be 
self-evident, we state the principle in the following formulation more suited for the issues under 
consideration: 

Time Invariance Principle (TIP): A physical law has the same mathematical form to 
every independent choice of observation time. 

This principle has an important implication to dynamical systems as laws of physical processes. 
To be precise, let x be the set of state variables and p be the set of parameters of a physical process. 
As a dynamical system, x changes in time t. Suppose an observation is made at t = and the state 
is xq. Another observation is made at time t > and the state is Xt. Then, as a physical law, Xt is 
governed by a function, denoted by Xt = 0t(xo, p), depending on the observation time t, the initial 
state xo, and the system parameter j>. As a default requirement, it must satisfy the unitary condition 



that is, with time increment 0, the law 0o leaves every state fixed. Now by the Time Invariance 
Principle, if another observation is made s > unit time later, the same function form {xt)s = 
'Psixt,p) must hold. Most importantly, the function 0f must satisfy the following group property 
and the unitary condition 



which together is referred to being TIP -conforming. That is, if an observation is made t time 
after the initial observation, and another is made s time later, then the result must be the same if 
only one observation is made s + t time after the initial observation. More generally, the state at 
s + t after an initial xq is the same state at s after an intermediate state Xt which is the state at t 
after the same initial a;o. A violation of this property that 0s+t(a^O;P) 7^ 4>s{4>t{xo,p),p) implies 
that either such an "experiment" is not reproducible, i.e., independent observation times inevitably 
lead to irreconcilable conclusions, or such a functional form does not govern the laws that the 
experiment is about to establish. 

An immediate consequence to the Time Invariance Principle is the following result. 

Lemma 1. If TIP -conforming dynamical system (j)t{x,p) is continuously differentiable at t = 
and any x in its domain of definition, then x(t) = (j)t{xQ,p) must be the unique solution to an 
initial value problem of a differential equation: 



4>o{xo,p) = Xo, 



{xt)s = (ps{xt,p) = (f)s{M^o,p),p) = (ps+t{xo,p) = Xs+t, and 0o(a;o,p) = xq, 



(1) 



dx(t) 
dt 



F^{x(t),p), x(0) = Xo, 
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where 



h=0 



is called the generating vector field ofcpf 

Proof. Because is differentiable and is TIP-conforming Qi, we have the following derivative 

dxjt) ^ 0t+fe(xo,p) - 0t(a^o,p) ^ ^.^ 0fe(0t(a^o,p),p) - MM^o^P)) 
dt h^o h h^o h 

= ^{(l>tixo,p),p) = F^ix{t),p), 
on h=o 

showing x{t) is a solution of the equation. Since p) is continuous differentiable in x because 
0j (x, p) is, the solution to the initial value problem is unique. □ 

We now conclude that Logistic Map does not model any population dynamics subject to time 
independent observations. More precisely, we have the following result. 

Theorem 1. There does not exist a continuously differentiable, TIP-conforming, 1 -dimensional 
dynamical system (pt{x, r) so that (pto {x, r) = Q{x, r) at any time to and for all xfrom any interval 
containing [0,1] for which Q is the Logistic Map and r is the intrinsic growth rate with r > 3. 

Proof By the preceding lemma, x{t) = (pti^Q^r) is the solution of an autonomous differential 
equation x' = F(f,{x) generated by 0. Since the system is 1 -dimensional, it does not allow periodic 
solutions. However, the Logistic Map has a period-2 orbit for r > 3 which would correspond to a 
periodic solution to the TIP-conforming flow if it were true that (x, r) = Q{x, r) for some to- 
A contradiction. □ 

This conclusion not only applies to the Logistic Map, but also to most other discrete maps in 
ecology. Table [H lists some popular discrete models in ecology. To be more precise, the same 
argument can be used to show the following. The generalized Beverton-Holt map is not TIP- 
conforming for 7 > 1 and large h. The same for the Bernoulli model for a > 1, the Richard 
map, the Ricker map for large r. Applying the same argument for 2-dimensional TIP-conforming 
functionals shows they are solutions to 2-dimensional autonomous differential equations which do 
not allow orientation reversing periodic orbits which occur to the Nicholson-Bailey map. 

For 3-dimensional or higher systems, the argument above for lower dimensional systems do 
not apply. However, here is a diagnostic test for possible TIP-nonconformity. More specifi- 
cally, we certainly assume that all biological processes are governed by physical laws that are 
TIP-conforming, allowing time-independent observation and verification on their states. Assume 
observation is made every unit of time and x„ is the state at time t = n. Because the state is 
TIP-conforming, we must have 

Xn = 0i(x„-i,r) = 0i(0i(x„_2, r),r) = ■ ■ ■ = 0"(xo,r) 

where the exponent stands for iterative composition. Therefore, 

(l)n{xo,r) = 0^(xo,r), 
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Table 1 : TIP-nonconforming Maps 



Generalized Beverton-Holt ([[30l[T9ll20]|) 

Bernoulli 
Logistic 
Richard (Il36l) 
Ricker (||37l) 



Nicholson-Bailey ( 113311 ) 



LPA (ttH) 



Leslie ( l|25l ) 



bNt 



r, 7 7^1 



*+l l+{hNt)i ■ 

Nt+i = aNtimod 1) 

Nt+i = Nt[l + r{l-Nt/K)] 

iVi+i = Nt[l + r(l - (Nt/Kr^] , m ^ 1 

Nt+i = Ntexp{r{l-Nt/K)) 

Nt+^ = Ntexp{-aPt) 
Pt+i = Ntil - exp(-aPt)) 



Lt+i 
At+i 



t+i 



= Lt{l- 
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that is, the nth iterative composition of 0i must have the same functional form as itself. This 
property can be used as a diagnostic test for probable TIP-nonconformity. For example, the nth 
iterate of the Logistic Map is a 2"-degree polynomial with evolving coefficients for each n. This 
implies that the map is very unlikely to be TIP-conforming because of the ever-changing functional 
forms of its iterates or its TIP-conforming functional would be extremely complex, in which case 
it is unlikely that such a complex functional happens to satisfy a stringent condition such as the 
TIP-conformity and at the same time arises from a relatively simplistic modeling exercise that 
is typical of most discrete modeling. This diagnostic test can be used to cast a serious doubt 
on the TIP-conformity of the model under consideration. Unlike the 1- and 2-dimensional maps 
discussed above, for which the preceding theorem provides a definitive means to determine their 
TIP-nonconformity, we can only conjecture based on the preliminary diagnostic test that the LPA 
map and the Leslie matrix are very unlikely to be TIP-conforming. The same can be said for 
all nonlinear models in cell-automata in games of life that they are very unlikely to model any 
physical processes subject to TIP-conformity. Without TIP-conformity, independent observations 
cannot verify nor establish such maps as models, theories, or laws. 

3. Discussions. The suggestion that biological research based on discrete maps is build on a shaky 
scientific ground inevitably leads to a few questions: One, is the Time Invariance Principle con- 
sistent with other known physical principles? Two, what are the TIP-conforming alternatives to 
discrete modeling in ecology? Three, can TIP-nonconforming maps be justified and under what 
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circumstances? Four, contradicting to predictions by all discrete, chaotic maps in theoretical ecol- 
ogy, why chaos is rare in the wild? We will exam these issues for the remainder of the discussion. 

Consistency With The Principles of Relativity. The Copernican idea that physical laws must be 
universal in space and time has guided many great theories in physics. Einstein's theories of special 
and general relativity are two of the most celebrated examples. The Time Invariance Principle is 
simply a corollary of the same idea that governs dynamical processes. 

One consequence of Einstein's theory of special relativity is that there is no absolute time. 
TIP captures this time-relativity aspect of his theory for the convenience of our discussion. As an 
example to make the point, consider two inertial frames with Frame 2 moving at a constant velocity 

V with respect to Frame 1 . Assume at the origin of Frame 1 there is an on-going dynamical process. 
Let ft{x) be the law deducted by observers of Frame 1 over a time interval t with initial state x. 
Due to the time dilation effect of special relativity for Frame 2, observers in Frame 2 will not see 
the same output ft{x) even though or precisely because they use the same clock time interval t. 
Instead, they will see /^(a;) for r = t/7 > t, 7 = a/1 — v'^/c^ with c being the speed of light. That 
is, both will see the same law but at two different observation times because of their best intention 
to use synchronized clocks at rest. This effect is equivalent to two observers in the same inertial 
frame using independent sampling times. 

At the center of the special relativity lies the Lorentz Transformation, relating space-time co- 
ordinates between two inertial frames. Let X = {t, x, y, z) be the space-time coordinate of Frame 
1 and Xy = Xy^y^^Zv) be the space-time coordinate of Frame 2 moving at a constant velocity 

V with respect to Frame 1, say along the same x-axis. Let X^ = L^{X) be the Lorentz Transfor- 
mation between the two coordinates. It is well-known that it satisfies this self-consistent Composi- 
tional Invariance Property that L„(X„) = Lu{Li,{X)) = L^{X) where w = {u + v)/{l + uv/c^) 
with c being the speed of light. That is, a third frame. Frame 3, moving at a velocity u relative 
to Frame 2 along the x-axis is a frame moving at a velocity w relative to Frame 1. The oper- 
ation (u, v) — > w defines the so-called Lorentz group. The Lorentz Transformation is one of 
the most well-known nontrivial and linear maps in physics that is compositionally invariant. The 
TIP-conforming group property is just a simpler kind of this more generalized Compositional In- 
variance Property. Stochastic matrixes (those which have non-negative entries and whose columns 
each sums to 1) form another well-known class of compositionally invariant linear maps. 

Most Dynamical Systems Should Be Modeled By Differential Equations. It is a well-known 
simple fact that if ipit, xq, p) is the solution of a differential equation x' = F{x, p) with t the time, 
xo the initial state at t = 0, and p the parameter, then it always satisfies the TIP-conforming group 
property 

for continuously differentiable right-hand side F. Continuous-time and probabilistic processes 
can be modeled by stochastic differential equations which also satisfy the TIP-conforming group 
property ([T]), c.f. Il26ll . In additions, the same group property is satisfied for hyperbolic PDF for 
age- structured populations, for parabolic PDF for reaction-diffusion and traveling wave phenom- 
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ena, and for delayed differential equations. In such cases, the states lie in some functional spaces 
which are infinitely dimensional. (Hence, by differential equations for the remainder of the dis- 
cussion we also mean to include such infinite dimensional equations with or without stochasticity, 
and as an extension by deterministic dynamical systems we also mean to include probabilistic 
processes modeled by stochastic differential equations for which it is the statistics in the means, 
variances, distributions, etc. of some state variables that become TIP-conforming, evolving deter- 
ministically.) 

TIP-conforming processes are not restricted to deterministic continuous-time processes only. 
In fact, true discrete probabilistic processes can be TIP-conforming. For example, the process 
of coin tossing has the same probability distribution satisfying 0„ = (pi with 0„ representing 
the distribution at the nth tossing, i.e., the same probabilistic law at any iterate of the process. 
Also, stochastic processes modeled by Markov chains are in general TIP-conforming because their 
transition matrixes have the same functional form for all iterates. 

The conclusion is that all (sufficiently differentiable, continuous-time) TIP-conforming laws 
are governed by differential equations. Hence, it is advisable to model all biological systems 
whose states are subject to independent time observation by differential equations in order to avoid 
the TIP-nonconforming trap that discrete maps can easily fall into. For true discrete probabilistic 
processes not modeled by differential equations, one must check the models TIP-conformity, for 
which stochastic matrixes represent one class of TIP-conforming, discrete, probabilistic models. 

The practitioner-dependent subjectiveness and arbitrariness of picking time increments in dis- 
crete modeling cannot be more apparent than modeling systems of various time scales. For ex- 
ample, when modeling a system of two species, such as an algae-zooplankton system or a plant- 
herbivore system, of which one operates at a faster time scale and the other operates at a slower 
time scale, what time step a discrete modeler should choose? If she picks the fast time scale to be 
her discrete time increment, she will end up fixing the slower species population as a parameter 
rather than an evolving variable, and missing out the population's temporal booms and busts. If 
she picks the slow time scale, she will end up aggregating the dynamics of multiple generations of 
the faster species, missing out its temporal booms and busts as well. Whatever choice she makes, 
it is very likely that her choice will not be honored by any independent modeler. Worse still, if 
her model is not TIP-conforming, which is extremely likely, her model is doomed. On the other 
hand, differential equation modeling does not have this time-scale misalignment problem of dis- 
crete modeling, which is easily dealt with by singularly perturbed differential equations, which 
will be discussed further later. 

TIP-Conformal Model — Derivation by One-Life Rule. Any textbook derivation of the Logistic 
Map as an one-species population model seems logically sound, yet it cannot be substantiated by 
independent time observation. The inescapable conclusion is that TIP must have been violated in 
all derivations. Two alternative fixes are proposed below. 

The first proposed fix follows the standard derivation of the Logistic Map with modification 
to the functional form of the per-capita growth of the species. The Logistic Map is the result of 
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assuming a linear functional for the per-capita growth as follows 



b — mxr 



where b is the maximal per-capita growth rate, and m is the mortality rate due to intraspecific 
competition. All empirical data (c.f. [[34l ). collected independently for different systems, with 
uncoordinated time increments, point to a density-dependent decreasing per-capita growth. That 
is, the decreasing monotonicity is qualitatively TIP-conforming for one-species per-capita growth. 
Although the linear functional is qualitatively TIP-conforming, it must have failed quantitatively, 
in particular at high population density. For example, if the density is x„ = (9 + 6) / m, then the per- 
capita growth is 6 — m(9 + b)/m = —9, implying paradoxically that each individual dies 9 times 
during the given interval of time, possible only for mythological cat. Furthermore, the per-capita 
growth any kth generation into the future, — Xn)/xn, fails to be strictly decreasing in x„ at 

hight growth rate r for any k > 2. This implies that the per-capita growth measured at a unit time 
interval t = 1 decreases in the population density, but does not when measured at, say, two units 
of time interval with t = 2. In other words, independent experiments would give contradicting 
outcomes in the per-capita growth functional if the Logistic Map were right, a not-so-surprising 
paradoxical effect of the map's TIP-nonconformity. Such inconsistencies are not limited to the 
Logistic Maps alone. In fact, they plague all one-dimensional maps from Table [H most of which 
are failed attempts to correct the Logistic Map. 
Our proposed fix assumes instead 

Xri4-1 X'fi b TTlXfi 



Xfi 1 ~1~ TflXyi 

Like the linear growth functional, it also decreases in x„, qualitatively TIP-conforming. More 
importantly, the per-capita growth is bounded below from —1, and approaches —1 only as x„ — > 
oo. That is, individuals die, but each dies at most once in any fixed period of time — a self-evident 
but both fundamental and universal principle for all organisms which is referred to as the One-Life 
Rule. As a result, the model is 

{b + l)Xn rXn „ . . 

Xn+1 = — := :— := 5i(a;„, r,m), 

1 + mXn 1 + mXn 

where r = b + 1 > 1. This results in the Beverton-Holt model which was first used by Beverton 
and Holt in 1956 for fishery studies (^), but not derived from the One-Life Rule as we did here. 

It is simple to check that the Beverton-Holt model is TIP-conforming — the composition of 
two B maps is another B map: 

S2(x„,r,m) := Bl{xn,r,m) = 5i(x„,,r^m(l + r)). 

In fact, one can demonstrate in general, 

Bk{xn,r,m) := B^{xn,r,m) = Bi{xn,r^,m -), 
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whose definition then can be extended to any continuous time t: 

Xt = Bt{xo, r, m) := Bi{xo, r\ m- ^) = — . 

r - 1 1 + m^— j-xo 

That is, any choice of fixed time increments leads to the same functional form. More importantly, 
the discrete Beverton-Holt map is simply the time-1 Poincare map of the continuous counterpart, 
and in this sense it cannot be regarded as a true discrete map. 

Also, as a consistency check, one can easily show that the per-capita growth over any time 
interval t is always a decreasing function of the initial population xq because 

xt-xo fTT-^Xo 
Xq 1 + m^Xo 

and that the One-life Rule is always satisfied since 

Xt-XQ 

lim = — 1. 

The Ricker model is an improvement over the Logistic Map in that it does obey the One-Life 
Rule. But it is not TIP-conforming and produces the same pathological chaos prediction for one- 
species population. 

TIP-Conformal Model — Derivation by Mass Balance Law. The second derivation can be 
best argued in terms of stoichiometry. It recognizes that an organism is a package of elemental 
elements, obeying the law of mass conservation. For example, let us use carbon (C) as a basic 
unit to measure an individual organism's biomass for a one-species system. Arbitrarily fix a time 
increment, say t = 1 for definitiveness. Let x„ and Xn+i be the numbers of individuals for the 
current generation and the "next" generation respectively. Let be the amount of C available in 
the interval, i.e., a constant flux in C. Let c be the amount of C that is needed during the period 
for each individual which is to make to the next generation, i.e., the per-capita maintenance cost in 
C. Let a be the efficiency rate, which measures the proportionality of the new generation that each 
individual of the current generation gives rise to for each unit of resource in C. It is the per-capita 
growth-to-consumption ratio. Then — cxn+i is the amount available for transition to the next 
generation, and the product of — cxn+i, a, and x„ gives the next generation's population: 

Xn+l = {N - CXn+l) X a X Xn- 

Simplify to obtain 

Xn+i = — with r = Na, m = ac, 

1 + mxn 

the same Beverton-Holt model obtained above. 

The Beverton-Holt model for one-species population is qualitatively consistent with all em- 
pirical studies cited in [[MllSTlfTSl . In particular, any initial non-zero population converges to an 
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equilibrium: 

lim Xn = lim E„(xo,r, m) = lim 5i(xo, r„,m„) 



n— ►oo 



r xq r — 1 

= lim jj— j — = , 

n^oo 1 + rn^-—^XQ m 

for which r > 1 as a defauk assumption. In the context of the stoichiometry for which r = 
Na, m = ca, we see that the greater the nutrient influx A^, the greater the stable equilibrium. The 
same holds for smaller per-capita maintenance cost c as well. The model predicts that prosperity or 
efficiency or both promote stability, not chaos, consistent with experimental findings such as [31 J. 

TIP-Equivalence Example — The Logistic Equation. As an illustration, consider the continuous- 
time Beverton-Holt model 

Xt = Bt{xo,r,m) 



1 I — 1 

1 + m-^xo 

Identical to the discrete case, it is also straightforward to check the TIP-conforming group property 

Bs+t{xo, r, m) = Bs{Bt{xo, r, m), r, m). 
Thus, the generating differential equation to which xt is a solution is obtained as 

dXt d I X| 

d ( Inr \ 

= -Tr^h{xur,m)\h=Q = Inr - m -Xt Xt, 



the Logistic Differential Equation! The per-capital growth rate is linear which can be arbitrar- 
ily negative at high population density. This does not violate the One-Life Rule which is for the 
per-capital growth in a fixed time interval rather than the instantaneous rate. The analysis above 
reaffirms a view that the Logistic Differential Equation is a good population model. By our argu- 
ment it is because it is TIP-equivalent to the Beverton-Holt model. 

Justifications of TIP-nonconforming Maps in Theoretical Ecology. Continuous models of food 
chains of three species or more can exhibit chaotic dynamics, which has been known since the ear- 
lier days of the chaos theory, see for examples, [[22l[T7ll2Tll32ll39ll24| . One particularly effective 
method to establish the existence of chaos in such models is the method of singular perturbation, 
see for examples, [|7l[8l|9l[T0l[ITl. At the so-called singular limit, some Poincare return maps are 
one-dimensional, nonlinear, and chaotic. However, such maps are not obtained by fixed time inter- 
val samplings. Instead, they are event return maps. For example, such a map may be defined when 
one of the predators reaches local maximums in population density, for which to occur the mo- 
ments in time cannot be independently chosen. Such event Poincare maps are different from fixed 
time-step Poincare maps in 3 critical ways: First, the sampling time for the former covers a con- 
tinuum range of interval, conditioned on the occurrence of the event. Second, the event Poincare 
maps are not TIP-conforming, but the fixed time-step Poincare maps are like their continuous time 
flows. Third, the fixed time-step Poincare maps are always at least 1 -dimensional higher than their 
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event counterparts. One-dimensional event return maps are of many types, including the unimodal 
type of which the Logistic Map is prototypical as demonstrated in the cited references. In conclu- 
sion, discrete event Poincare maps are used as auxiliary means to understand continuous models 
which generate them. They alone do not model the underlining process in any time independent 
fashion. 

Discrete maps also play an irreplaceable role in numerical approximations. In fact, the Lo- 
gistic Map is a discretization of the Logistic Differential Equation. It is not TIP-conforming for 
any r but it is a good approximation of the TIP-conforming continuous model when r is near 1 . 
Similarly, TIP-nonconforming maps generated as numerical schemes to approximate continuous 
models do serve useful and important purposes in theoretical and practical applications. But they 
are relevant only within their realistic ranges. For example, it has always been an unjustifiable 
extrapolation to large r of the Logistic Map that it becomes artificial and problematic. More ad 
hoc still than the Logistic Map, the Ricker Map is not a discretization of any known differential 
equation, further removed from TIP-conformity. Similarly, a Leslie matrix should have been de- 
rived as a discretization of the linearization of its PDE counterpart of the age- structured population 
intended by the matrix, which unfortunately was not always the case. Such a discretization of the 
linear PDE should impose constraints on the time step and the age increment. Outside such con- 
straints, the Leslie model may become problematic as is the case for the Logistic Map. With the 
above rather thorough analysis of the Logistic Map, one implication seems hard to miss that dis- 
crete models without underlining TIP-conforming origins only have limited if not all questionable 
scientific values. 

Enrichment and Efficiency Stabilization Principles. As pointed out in the introduction that em- 
pirical studies do not support the hypothesis of chaotic one-species population at high reproductive 
per-capita rate, or high efficiency rate. In fact, the Beverton-Holt model implies diametrically the 
opposite. For higher dimensional systems, the same opposing dichotomy existed, and it can also 
be reconciled by TIP-conforming models. 

On one hand, chaotic dynamics do occur in models for three species or more in food chains 
and webs (JH Il7l|2ll|32ll2l[35lll|9l[l0lllllll|5l Il3l), as well as in laboratory models of 
3-dimensional systems (f6]|). Also, almost all models exhibit the same paradoxical effect that 
high reproductive efficiency leads to chaotic population dynamics. Recent analysis from [[T2l[T3l 
concluded that this "Chaos Paradox" is an artifact of the Malthusian exponential growth model 
for populations. In addition to its projected unbounded growth fallacy, the Malthusian model, 
P' = rP, also violates the One-Life Rule: {P{t) — Pq)/ Pq = e*"* — 1, not going to —1 as Pq oo. 
As pointed out in [fT2l[T3l . this Malthusian pathology hides in almost all continuous models in the 
literature, and paradoxical results are inevitable. The Enrichment Paradox ( ifBSl ). the Biological 
Control Paradox dllSl ). the Competition Exclusion Principle ([HI), and the Chaos Paradox are of 
the most notorious. 

On the other hand, however, assuming logistic growth for all species provides a sufficient rem- 
edy to all these paradoxes, the Chaos Paradox in particular. For example, consider the following 
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food chain equations 

X =X {h -d^- m,X) - ""f Y 

1 + hiaiX 

Y =Y { -d,- m,Y\ - Z (2) 

f ha2Y 

Z =Z ; — — d'i — rri'iZ 

\l + h2a2Y 

The result of [fT3l shows that all chaotic attractors bifurcate into either limit cycles or steady states 
as the top-predator's reproductive efficiency parameter 63 increases, and that the limit cycles further 
bifurcate into steady states if the predator's reproductive efficiency parameter 62 also increases. In 
fact, one can show that all chaos attractors must bifurcate into a steady state by increasing any two 
of the three reproductive efficiency parameters in 61, 62, &3- In other words, efficiency promotes 
ecological stability, a result of evolution if the hypothesis holds that evolution promotes species 
efficiency as a survival fitness. 

TIP-conformity is a necessary but not a sufficient requirement for physical laws. Comparing 
to most studied in the literature, the food chain model ^ above is significantly better because 
of two features incorporated into the model. First, with the inclusion of parameters mi, m2, m^, 
intraspecific competitions are taken into consideration for all species which leads to the logistic 
growth rate for individual species ( flOl ) which in turn conforms to the One-Life Rule. Second, 
the Rolling Type II predation functional form ( ll23l ) is used for all predators. The importance of 
this particular form lies in the fact that it is TIP-conforming from its own derivation, and more 
generally it satisfies the Composition Invariance Principle. To see this, we recall from Rolling's 
original mechanistic derivation that r^, = aX is the number of prey, X, encountered in a unit 
time by one predator, Y, or the encounter rate, and h is the handling time per kill of the prey, 
and r = -j^^ = i^^^x number of handled kills in one unit time by one predator, or the 

predation rate. This functional form is TIP-conforming for the following reasons. If one breaks 
the handling time down to, say, per-killing time hk, and per-consuming time he with h = hk + he, 
then Rolling's derivation will give rise to the kill rate = j:^^^^ as a function of killing time and 
encounter rate, and the consumption rate Tc = ^^^^/^ as a function of consuming time and kill rate. 
It is straightforward to check that 

_ fk _ l+/ifcre _ "^e 



" l + Vfc l + h,j^^ l + {hk + K)r; 

That is, the rate function has the same mathematical form regardless the temporal cut-off or defi- 
nition of handling time, the essence of TIP-conformity and compositional invariance. 

In conclusion, theoretical and qualitative predictions of the food chain model ^ are consis- 
tent with those of the Beverton-Rolt model as well as relevant empirical findings that enrichment 
and efficiency promote ecological stability. The consistency may not be coincidental, because all 
compartmental constituents of the model are mechanistically TIP-conforming. 
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Concluding Remarks. Empirical data almost always are collected at discrete times. Discrete 
modeling is an intuitive response to that reality to fit discrete data by discrete models. However, a 
discrete model has little to say about data collected at different discrete times of the same process if 
the process permits. If it is not TIP-conforming, it does not model the underlying process subject to 
time independent observation. This may underlie many attempts via stochastic inclusion to discrete 
modeling, attributing noise or stochasticity as the chief cause of the irreconcilability between a 
theory and reality when in fact TIP-nonconformity of the theory may have been the problem. 
TIP-nonconforming event maps are secondary structures of inherently higher dimensional TIP- 
conforming differential equations. They rarely have a closed-form formula with system parameters 
in plain sight for meaningful manipulations because of the aggregating procedures that produce 
them. Even in such cases, they are not closer in capturing the underlying physical laws for the 
processes than the TIP-conforming differential equations that model the processes. Given all these 
considerations, this paper advocates a typical approach to use continuous models to fit discrete 
data. Such models are open to the scrutiny of all observations carried out at any discrete times. 
This approach makes sure the models are necessarily consistent in its internal and conceptual 
construct, allowing the modelers to modify and to refine them within the realm of TIP-conformity. 

Our TIP-equivalence result for fixed time-step Poincare maps and differential equations implies 
that 1 and 2 dimensional TIP-conforming maps cannot be chaotic because 1 and 2 dimensional 
differential equations of continuously differentiable vector fields cannot be chaotic. Equivalently, 
chaotic 1 and 2 dimensional maps must be TIP-nonconforming, and at the best arise as event 
Poincare maps of 3 or higher dimensional differential equations. As a result such maps do not 
model any physical processes at a time-independent fashion. Hence an ecological conclusion can 
be made unequivocally that single- and two-species population dynamics cannot be chaotic. That 
the controversy has lasted this long was due to the combination of a few understandable factors. 
To name a few obvious: First, the derivations of all popular discrete ecological models seemed 
logical, but TIP-nonconforming nevertheless. Secondly, because of their TIP-nonconformity, all 
predictions could not be independently and objectively reproduced, leading to the inevitable con- 
fusing state between a seeming reasonable theory and an uncompromising reality. Thirdly, the field 
irreproducibility of all low dimensional chaos theory was conveniently masked by the inherited 
unpredictability of all chaotic systems. And fourthly, the irreconcilability was also conveniently 
masked by a noisy reality that is for most biological experiments and observations. 

Comparing to differential equations, discrete maps are easier to teach, easier to do research 
with. But we should not compromise the Time Invariance Principle just for the simplicity appeal 
of discrete modeling. TIP-conformity is the minimal necessary condition a conceptually consis- 
tent model must satisfy. More importantly, the requirement is fundamental to all branches of 
science, governing the reproducibility of experiments. Because of these reasons, usage of TIP- 
nonconforming maps is difficult to justify in most circumstances. This conclusion has some 
important implications to both research and training: Both past and future researches based on 
discrete models must be scrutinized against their TIP-conformity and be justified for their TIP- 
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nonconformity. The subject of discrete modeling may have to be de-emphasized in the classrooms 
and be viewed through the lens of TIP-conformity. On the other hand, training in calculus and 
differential equations must be further enhanced and greatly emphasized for future generations of 
theoretical biologists. 
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